Statistics and noise in a quantum measurement process 
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The quantum measurement process by a single-electron 
transistor or a quantum point contact coupled to a quantum 
bit is studied. We find a unified description of the statistics of 
the monitored quantity, the current, in the regime of strong 
measurement and expect this description to apply for a wide 
class of quantum measurements. We derive the probability 
distributions for the current and charge in different stages of 
the process. In the parameter regime of the strong measure- 
ment the current develops a telegraph-noise behavior which 
can be detected in the noise spectrum. 

Introduction. The long-standing interest in funda- 
mental questions of the quantum measurement received 
new impetus by the experimental progress in mesoscopic 
physics and increasing activities in quantum state engi- 
neering. The basic idea is to use as meter a device, able 
to carry a macroscopic current, which is coupled to the 
quantum system such that the conductance depends on 
the quantum state. By monitoring the current one per- 
forms a quantum measurement, which, in turn, causes 
a dephasing of the quantum system |^-|^ . The dephas- 
ing has been demonstrated in an experiment of Buks et 
al. Q where a quantum dot is embedded in one arm of 
a 'which-path' interferometer. The current through a 
quantum point contact (QPC) in close proximity to the 
dot suppresses the interference. However, since passing 
electrons interact with the current only for a short dwell- 
time in the dot, the meter fails to distinguish between 
two possible paths of individual electron; only a reduc- 
tion of interference has been observed. This situation is 
referred to as a weak measurement. 

For a strong quantum measurement a different setting 
is needed, where a closed quantum system is observed by 
a meter. Then a sufficiently long observation may provide 
information about the quantum state. This situation is 
realized when a single-electron transistor (SET) is cou- 
pled to a Josephson junction single-charge box, which for 
suitable parameters serves as a quantum bit (qubit) [^,^. 
The analysis of the time evolution of the density matrix 
of the coupled system demonstrates the mutual influence 
between qubit and meter, i.e. measurement and dephas- 
ing §. 

The measurement process is characterized by three 
time scales. On the shortest, the dephasing time r^p, 
the phase coherence between two eigenstates of the qubit 
is lost, while their occupation probabilities remain un- 
changed. Later measurement-induced transitions mix the 
eigenstates, changing their occupation probabilities on 
a time scale t^ix > and erasing information about 



the initial state of the qubit. The origin of the mix- 
ing is that the charge operator (the measured quantity) 
and the qubit's Hamiltonian do not commute. The third 
time scale appears in the dynamics of the current in the 
SET. Consider the probability distribution P(m, t) that 
m electrons have tunneled through the SET by time t. It 
was shown that after a certain time tmcas information 
about the state of the qubit can be extracted by reading 
out m. As expected from the basic principles of quantum 
mechanics the observation of the qubit first of all disturbs 
its state. Hence imcas > t^- The measurement process is 
only effective if the mixing is slow, i„iix ^ ^moas- In the 
opposite limit of strong mixing the information about the 
qubit's state is lost before a read-out is achieved. 

The distribution function P{m, t) describes the statis- 
tics of the charge which has tunneled. The distribution 
of the current in the SET p(I, t) and current-current cor- 
relations require, furthermore, the knowledge of corre- 
lations of the values of m at different times. In earlier 
papers on the statistics in a SET or a QPC the 
behavior of P{m, t) at times shorter than t^i^ was de- 
rived. Effects due to the additional knowledge acquired 
by an observer ||l^, and the possible influence of the 
wave-function collapse on the monitored current [|j were 
also discussed. Here we develop a systematic approach, 
based on the time evolution (Schrodinger equation) of 
the density matrix of the coupled system. This approach 
allows us to study averages and correlators of the current 
and charge. Since, due to shot noise, instantaneous val- 
ues of the current fluctuate strongly, we study the current 
/, averaged over a finite time interval At. We calculate 
the mixing time in a SET and derive analytic expressions 
for p{I, At, t), as well as P(m, t), vahd on both short and 
long time scales. We study the noise spectrum of the 
current and find that in the limit of strong measurement 
(^mcas < ^mix) the long-time dynamics is characterized by 
telegraph noise, with jumps between two possible current 
values, corresponding to two qubit's eigenstates. 

The results are of immediate experimental interest. A 
recent experiment demonstrated the quantum coherence 
in a macroscopic superconducting electron box |ll| , but 
the coherence time was limited by the measuring device. 
The SET-based measurement should extend the coher- 
ence time, which combined with experimental progress 
in fast measurement techniques [ p^ should increase the 
maximum number of coherent quantum manipulations. 

Master equation for the measurement by a SET. The 
system of a qubit coupled to a SET is shown in Fig. [|. 
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The qubit is a Josephson junction in the Coulomb block- 
ade regime. Its dynamics is limited to a two-dimensional 
Hilbert space spanned by two charge states, with n — 
or 1 extra Cooper pair on a superconducting island. The 
island is coupled capacitively to the middle island of the 
SET, influencing the transport current. The SET is kept 
in the off-state during manipulations on the qubit 0, 
with no dissipative current and no additional decoher- 
ence. To perform the measurement, the transport volt- 
age is switched to a sufficiently high value, so that the 
current starts to flow in the SET. 
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FIG. 1. The circuit of a qubit and a SET used as a meter. 



The Hamiltonian of the system is given by |lj 



Ti- ~ HsET + + Ti-T + 'Hip 



(1) 



The first term is the charging energy of the transis- 
tor, quadratic in the charge A^e on the middle island, 
■^SET = -EsetI^ — ^o)^- The induced charge N^e is 
determined by the gate voltage Vg and other voltages 
in the circuit. The Hamiltonian of the qubit, = 
"Hqh + NTiint, includes the Hamiltonian of the uncoupled 
qubit Hqb and the Coulomb interaction with the SET 
NHint (split into two terms for later convenience.) In the 
basis of the qubit 's charge states they are given by Hqh — 
^{Eqh^z - E,ja^) and NHint = EintNa^- The charging 
energy scales i?sETj ^-qb and Ei^t are determined by ca- 
pacitances in the circuit, while Ej is the Josephson cou- 
pling. We consider the eigenstates of Tiqb, |0) and |1), as 
the logical states of the qubit. In this basis, Hqh is diag- 
onal, with the level spacing AE = {E'j -f -Bqb)^^^, while 

* where e}^, = E,,,tEqh/AE and 
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E^^^ = EintEj/ AE. The term TL^p describes the Fermions 
in the island and electrodes of the SET, while TLt gov- 
erns the tunneling in the SET. Here we assume weak 
coupling to the environment, with relaxation slower than 
the SET-induccd mixing. The opposite limit is dicussed 
in Ref. |f|. 

The full density matrix can be reduced to Pj%i™i (t) by 
tracing over microscopic degrees of freedom and keeping 
track only of N and m, the number of electrons which 
have tunneled through the SET. Here i,j = 0, 1 refer to 
a qubit's basis. A closed set of equations can be derived 
for p'^{m) = p'j^^, the diagonal entries of the density 
matrix in N and m p^ . Solving these equations, we 
analyze the evolution of the reduced density matrix of 
the qubit, Q'^^{t) = YliN mP^Ni''^^'^)^ "^^ll ^ P{m,t) = 
J2n j p'n (™' ^^'^ other statistical characteristics of the 



current in the SET. 

At low temperatures and transport voltages only two 
charge states of the middle island of the SET, with N 
and N + 1 electrons, contribute to the dynamics. Trans- 
lational invariance in m-space suggests the Fourier trans- 
formation p^^{k) = ^~''^™P7v("^)- Expanding in the 
tunneling term to lowest order, we obtain the following 
master equation (cf. Refs. [Rjl^): 
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i ["^qb, Pn] 

i [^qb + Hint, PN+i] 



Here the operators 



f L/5 = TlP + TTUl [Hint: p] + 
^RP = ^RP - T^OIR [H.int,p\ + 



(2) 



(3) 
(4) 



represent the tunneling rates in the left and right junc- 
tions, with aL/R = ^k/ i&T^^ R]^/^) being the tunnel con- 
ductance of the junctions in units of the resistance quan- 
tum i?K = h/e^. The rates are fixed by the potentials 
PL, PR ^ PL + Vtr of the leads: F^ = 2TraL[pL - (1 - 
2No)Eset] and F^ = 2^afl[(f - 2iVo)SsET - Pr]- They 
define the tunneling rate F = ^L^R/i^L + ^r) through 
the SET. The last terms in Eqs. (^,^ make these rates 
sensitive to the qubit's state. 

The initial condition at the beginning of the measure- 
ment, written down in the logical basis. 



p%{m,t = 0) = 



|a|2 ab* 
a*b 



WN 



mO 



(5) 



describes the qubit in a pure state a|0) I) and the SET 
in the zero-voltage equilibrium state. One can assume 
that wn = 1 and wat+i = 0. 

Reduction of the master equation. In general the dy- 
namics of the qubit's density matrix g, described by the 
master equation (|^), is complicated since dephasing (de- 
cay of the off-diagonal entries) and mixing (relaxation 
of the diagonal to their stationary values) may occur on 
similar time scales, « imix- However, under suitable 
conditions the mixing is slow, which is the prerequisite 
for a measurement process. This is the case, if the qubit 
operates in the regime with dominant charging energy: 



l-E-qbl, \Eqh + 2ii'int| S> Ej 



(6) 



Weak {Eint ^ AE) or strong {Eint ~ AE) coupling to 
the SET can be considered. In the latter case a faster 
measurement is achieved (see Eq. (pi])). 

We first analyze the dynamics without mixing, i.e. we 
put E^^ — 0. In this case the time evolutions of for 
the four different pairs of indices ij are decoupled from 
each other, each being characterized by two eigenmodes. 
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The absence of mixing, further, imphes the conservation 
of occupations of the logical states g", for i = 0, 1, and we 
find two 'Goldstone' modes for fc ^ 1, with eigenvalues 



1 



A!^(fc) w ipfc- -/T'fc 



(7) 



Here T' = r^r*^/(r*^ + T^^) are the tunneling rates 
through the SET, T^^^ = TL±TTaLEl^ and r°/^ =TrT 
7rafli?Jjj^ are the tunneling rates in the junctions for two 
logical states (cf. Eqs. (||)), and P = 1 - 2rV(P^ +r}j) 
are the Fano factors responsible for the shot noise re- 
duction. The other two eigenmodes decay fast, with the 
rates A*:! w -(P^ + 

The analysis of the eigenvalues, A^^(fc) = [A?[?(— /c)]*, 
of the four off-diagonal modes in ij, reveals the dephasing 
time of the qubit by the measurement, i.e., the decay 
time of = p°j}{k = 0) + = 0)- I* gi^^^ 

by T-i = ATEl^/iTL + TRf if « Tl + Tr, and 
= wjvTl -I- wn+iTr in the opposite case. 

The picture is modified by the mixing at finite -Bj^jj. 
We find that the mixing may be treated perturbatively if 
I Aj?^ I ^ E-;^^^ , which turns to be the case if the condition 
holds. Then, in the second order, the degeneracy be- 
tween the long-living modes (0) is lifted and the long-time 
evolution of the occupations p"(fc) = p%{k) + p5v+i(^) 
given by a reduced master equation, 



d 
di 



P"(fc) 



Mrcd 



pO"(fc) 
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-0 



^^^--'^~[ AV(fc) 



For the mixing rate, Fmix, we obtain 

^ mix ~ 

AE^iAE + 2eI^)^ + [TrAE + T^AE + 2eI 



(8) 



(9) 



(10) 



To understand the role of the mixing we assume first 
Fmix = in Eqs. (|8|J^). Then, for the initial condition 
(I) we obtain p°°{k) = \a\^e^+^'''> *, p"(fc) = |6|2e^+ *, 
and P{k, t) = P{m, t) e''^"" = p°° (k) + p^\k) . From 
this we obtain the distribution P{m,t), which evolves 
from a peak S{m) at i = into two peaks with weights 
|ap and moving in m-space with velocities and 
F , and with widths growing as ^2fYH. The peaks 
separate after a time 



FO -Fi 



(11) 



Thus measuring the charge m after tmcas constitutes a 
strong quantum measurement 0. However, at longer 
times t > Fjjj?^ the mixing spoils this picture. In par- 
ticular, the occupations of the logical states relax to the 



equal-weight distribution: (t) — g^^ (t) oc exp(— Fmixi). 
Therefore the two-peak structure appears only in the in- 
terval between tmeas < i < ^mL- ^^'^ measurement can 
be performed only if t^eas < T^^ly, 

The measurement takes longer than the dephasing, 
^moas 3> Ty . Such measurement can be called non- 
efficient |in]: the information about the qubit is con- 
tained in the SET already after r^, but can be read out 
from the current only later. 

The quantum measurement with a QPC can be de- 
scribed in a similar way. The Coulomb interaction of 
the qubit with the current results in two tunneling rates 
po/i — p -|- SY/2 for two qubit's states. Tracing out mi- 
croscopic degrees of freedom one arrives at the master 
equation ||^ for the Fourier transform of the density ma- 
trix p*-'(m), which can be rewritten as 



tp+^ST[a,,p]+ (e"= 
-(4T^)-ie"=[a„[a.,p]] 



1) 



(12) 



One can show that r,. 



i(Vro - %/F)2 is the decay 
rate of For p" the eigenvalues are given by (0), with- 
out Fano factors. The measurement time ( pT| ) and the 
dephasing time coincide, implying the 100% efficiency. 

Under the condition Ej <C max(Ai!^, r""'^) the per- 
turbative treatment produces the reduced master equa- 
tion (|^,^), with the mixing rate 



E 



1 + AE^ tI 



(13) 



A phenomenon, termed the Zeno (or watchdog) effect, 
can be seen in the Hmit ^ AE: the stronger 
is the measurement, quantified by t~^, the weaker is the 
rate of jumps between the eigenstates, Fmix ~ EjT^. 

The analysis of the SET mixing rate (|l^) in terms of 
the Zeno effect is more complicated. The rates t^^ and 
Fmix depend in this case on several parameters and no 
simple relation between Fmix and is found. However, 
in the regime i^int Tl ^ Tr <s:^ \'^E\, these two rates 
change in opposite directions as functions of F^, ^ F^, 
which is reminiscent of the Zeno physics. 

Statistics of charge and current. The results of this 
section apply to the SET and QPC alike. The statisti- 
cal quantities studied depend on the initial density ma- 
trix (||): e.g., P{m,t) = P{m,t \ po). In the two-mode 
approximation (HJ^) this reduces to a dependence on 
70 = - = |a|2 - We solve Eq. (§ to ob- 
tain P(r7i,i |7o) = Trqb[J7(m, i)po], where U{m,t) is the 
inverse Fourier transform of U{k,t) = exp [Mi.ed(^) t]- If 
po/i ^ f±<5F/2 are close, the resulting distribution is 

-<5m^/2/rt 

P(m,^|7o) = ^P(m-,5m,t|7o) . (14) 
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The first term in the convolution (jlj) contains two delta- 
peaks, corresponding to two qubit's logic states: 



'm — rt 1 
P(m,t\jo) = Ppi I , -r,„ixi 



70 



6n/2 ' 2 

+ e-^""-*'/^\a\^S{m-r"t) + \b\H{m-rh)] . (15) 

On the time scale imix = r„iix peaks disappear; in- 
stead a plateau arises. It is described by 

+ (1 + 7ox) /i (rv/l-x2) /v/r^} , (16) 

at |a;| < 1 and Ppi — for |x| > 1. Here Iq, Ii are the 
modified Bessel functions. At longer times the plateau 
transforms into a narrow peak centered around m = Tt. 

The Gaussian in Eq. ([l^ ) arises due to shot noise. Its 
effect is to smear out the distribution (see Fig. ^). 

We also calculate P2{m,t;m' ,t' \po), the joint proba- 
bility to have m electrons at t and m' electrons at t'. 
This allows us to obtain the probability distribution 



p(/, At,t|po) 



P2{m + IAt,t + At;m,t\po) (17) 



of the current / = /j*^'^* I{t')dt' averaged over the time 
interval At. The evolution is Markovian, and we obtain: 
P2(m, t;m + Am, t + At) = Trqb_[C/(Am, At) C/(m, t)po] 
for At > 0. The derivation of p{I, At, t) thus reduces to 
the calculation of the charge distribution (14) for differ- 
ent initial conditions: 



p(/,At,t|7o) = P(m = /At,At| 



*7o) . (18) 



The behavior of p(/. At, t) is shown in Fig. |^. A strong 
quantum measurement is achieved if tmcas < At < t^i^, 
at times t < tmix (see Fig. 0b). In this case the measured 
current is close to either r° or F^, with probabilities |ap 
and l&p, respectively. At longer times a typical current 
pattern is a telegraph signal jumping between F" and F^ 
on a time scale tmix- If At <^ tmcas the meter does not 
have enough time to extract the signal from the shot- 
noise background. At larger At the meter-induced mix- 
ing erases the information, partially {At ~ tmix, Fig. ^) 
or completely {At ^ tmix, Fig. ||d), before it is read out. 

The telegraph noise behavior is also seen in the current 
noise. Fourier transformation of the correlator 



{I{t)I{t')), 



E 



mm' dtdt' P2{Tn,t;m' ,t' \ po) (19) 



gives in the stationary case the noise spectrum. 



Si{uj) = 2eVr + 



as the sum of the shot- and telegraph-noise contributions. 
At low frequencies cu <C Fmix the latter becomes visible 
on top of the shot noise as we approach the regime of the 
strong measurement: S'toiograph/S'shot ~ 4tmix/tnioas- 



P{m/t,t)-t 




p(I,M,t) 




p{I,At,t) 




^meas ^ 



p{I,M,t) 




At>t, 



FIG. 2. The probability distributions of the charge (a) and 
current (b-d). P{m,t) in (a) is rescaled, so that the peaks do 
not move. The time-axis scale is logarithmic. 

To conclude, we have developed the master equation 
approach to study the statistics of currents in a SET or 
a QPC as a quantum meter. We evaluate the probability 
distributions and the noise spectrum of the current. 

We acknowledge discussions with Y. Blanter, 
L. Dreher, S. Gurvitz, D. Ivanov, and A. Korotkov. 
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